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principally  the  solar  phase  angle,  into  the  prediction  of  aircraft  luminance 
The  test  consisted  of  a comparison  of  model  prediction  with  field  test  data 
consisting  of  photographic  imagery  of  a low-altitude  aircraft. 

The  photographic  data  were  reduced  through  two-dimensional  scanning 
densitometry  to  yield  mean  silhouette  luminance.  These  values  were 
compared  to  values  predicted  by  a diffuse  sphere  model  with  the  comparison 
including  a full  account  of  atmospherics.  A discussion  of  the  comparison 
is  included. 


*J&L IMJ*. 

Wlte  VicW 

1 2 » >“'*■  c 

1 noc  c 

\ USHW38KO 

1 j'dsnfiMTio* 

i :: 

1 t0>-’ 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS%PAGEfltTien  0»f<*  Fnlerf ./) 

■ ill I H 


SUMMARY 


Visual  search  for  approaching  aircraft  by  ground  observers  is  a part 
of  many  system  performance  evaluations.  As  a result,  several  mathematical 
models  exist  that  quantify  the  observers'  ability,  usually  in  a probabilistic 
fashion.  Basic  to  all  the  models  are  certain  fundamental  measures  such  as 
mean  luminance  of  the  aircraft  silhouette  as  portrayed  in  aspect,  background 
luminance  against  which  it  is  presented,  and  angular  size  apparent  to  the 
searching  observer. 

Presented  in  this  report  are  the  results  of  a study  to  test  one 
particular  model  for  aircraft  luminance.  Specifically,  the  model  is  simply 
to  represent  the  aircraft  as  a diffuse  (Lambertian)  reflecting  sphere.  This 
allows  a sun  angle  dependence  of  sunlit  and  shadow  regions  that  contribute 
to  total  luminance  to  be  included  as  a first  approximation.  This  representa- 
tion was  combined  with  surround  luminance/atmospheric  model  to  yield  a total 
model  for  comparison  to  available  photographic  data  collected  in  a field  test. 
The  surround  luminance  model  includes  the  ground  and  sky  as  sphere  illumi- 
nators, as  well  as  the  direct  solar  component.  The  atmospherics  account  for 
reduction  of  illumination  according  to  solar  slant  path  through  the  total 
air  mass  and  transport  between  the  aircraft  and  camera  that  yielded  the  data. 
Although  atmospherics  were  not  regarded  as  very  important  in  the  field  test, 
as  much  detail  as  possible  was  included  in  the  analysis  to  assure  it  would 
not  be  an  issue. 

The  data  for  model  comparison  consisted  of  16  mm  imagery  of  a low- 
altitude  jet  aircraft  collected  with  a camera  system  corrected  to  a photopic 
response.  The  imagery,  which  was  collected  over  all  azimuths,  was  supported 
by  radar  coordinates  of  the  aircraft  such  that  required  geometry  would  be 
available.  Also  included  in  the  imaging  were  luminance  control  panels, 
the  luminance  values  of  which  were  measured  simultaneously  with  film  exposure, 
thereby  providing  absolute  film  calibration. 

The  imagery  was  analyzed  with  a scanning  two-dimensional  densitometer 
with  16  levels  approximately  proportional  to  film  transmission.  Thus,  for 


iii 


* 


each  frame  analyzed,  frame  areas  could  be  measured  for  each  level,  which, 
in  conjunction  with  the  calibration  data,  allowed  mean  aircraft  luminance 
(apparent  at  the  camera)  to  be  computed. 

The  data/model  comparison  yielded  significant  departure  of  the  data 
from  the  model  at  certain  solar  phase  angles.  These  data  points  all  possessed 
sun  angle  geometry  together  in  the  aircraft  surface  normal  directions  consis- 
tent with  solar  glint.  Other  than  these  data  points,  a trend  of  decreasing 
luminance  with  phase  angle  was  observed,  but  weaker  than  model  prediction. 
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Section  1 
INTRODUCTION 

The  prediction  of  ranges  at  which  aircraft  can  be  detected  visually 
by  ground  and  airborne  observers  is  of  perennial  interest.  To  make  such 
predictions  with  any  degree  of  confidence,  information  regarding  the  intrinsic 
contrast  of  the  aircraft  against  its  background  is  required.  Because  of  the 

i 

multiplicity  of  types  of  aircraft,  of  background  against  which  they  might 
be  viewed,  of  possible  lighting  distributions,  and  of  the  difficulty  of 
measurement  of  variation  in  contrast  of  a target  during  flight,  very  little 
systematic  reliable  contrast  data  exists. 

This  report  presents  the  results  of  a study  to  test  the  validity 
or  appropriateness  of  one  particular  mathematical  model  for  aircraft  luminance. 
Visual  detection  considerations  (e.g.,  the  role  of  resulting  contrasts  in 
detectability)  were  of  no  concern  in  the  study.  Specifically,  the  approach 
is  one  of  comparing  model  prediction  to  available  data  to  reach  a conclusion 
on  model  utility. 

The  model  tested  is  the  representation  of  an  aircraft  as  a diffuse 
• (Lambertian)  reflecting  sphere  for  the  purpose  of  computing  mean  aircraft 

luminance  over  the  silhouette.  Although  real  aircraft  are  neither  diffuse 
nor  spherical,  such  a representation  does  allow  some  account  of  the  variability 
with  sun  angle  of  the  contributions  due  to  sunlit  and  shadowed  regions  with- 
out undue  model  complexity.  Calspan  has  in  the  past  used  such  a model  as  an 
indicator  of  contrast  variability  to  assist  in  selection  of  contrast  values 
appropriate  as  inputs  to  a visual  detection  model.  Others  have  also  used  such 
a model,  in  particular,  Nagel  proposes  its  use  to  represent  "full  form" 
objects  generally,  and  his  paper  contains  considerable  detail. 


1.  M.R.  Nagel,  "A  Model  for  the  Inherent  Contrast  Conditions  in  Full-Form 
Objects,"  Deutsche  Forschungs-  und  Versuchsanstalt  fur  Luft-  und 
Raumfahrt,  D-8031  Oberpfaffenhofen,  Germany. 
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Two  sets  of  data  were  considered  for  model  comparison:  (1)  some 

results  from  Project  COMPASS  GHOST  in  the  form  of  ground-to-air  photography 
of  low-altitude  jet  aircraft,  and  (2)  imagery  and  support  data  collected 
during  a flight  test  by  MICOM/USAADS  within  this  study  itself.  An  attempt 
was  made  to  analyze  the  COMPASS  GHOST  imagery,  but  this  was  abandoned  due  to 
insufficient  absolute  film  calibration. 

Included  in  the  sections  that  follow  are  a description  of  the 
diffuse  sphere  model,  a discussion  of  the  flight  test,  and  a mathematical 
account  of  the  effects  of  atmospherics  also  needed  when  the  comparison  is 
made.  Finally,  the  model/data  comparison  is  made,  and  conclusions  are 
. drawn . 
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Section  2 

DIFFUSE  SPHERE  MODEL 


The  details  of  the  application  of  a diffuse  reflecting  sphere  for 
predicting  aircraft  brightness  are  given  in  this  section.  The  discussion  is 
limited  to  the  local  reflection  process;  viz.,  the  relationship  between  flux 
incident  upon  the  sphere  and  the  reflected  flux.  Other  aspects  necessary  to 
ready  the  model  for  comparison  to  field  test  data,  such  as  brightness 
propagation  through  the  atmosphere,  are  included  in  later  sections  where 
appropriate.  Finally,  it  is  important  to  understand  that  the  use  of  the 
diffuse  sphere  here  is  not  the  same  as  saying  that,  from  the  viewpoint  of 
visual  detectability,  the  aircraft  is  modeled  as  a sphere.  It  is  used  only 
to  establish  average  aircraft  brightness  with  a first-order  dependence  on 
sun  angle;  other  detection  parameters,  such  as  aircraft  size,  can  be  obtained 
from  more  appropriate  sources. 

The  average  brightness  of  a diffuse  reflecting  sphere,  as  it 

depends  on  geometry  assuming  remote  point  illumination,  has  been  computed 
2 

by  Russell.  We  can  write: 

B = ^ [5  in  oc  ■+■  (rr  - oc)  cos  <aTj 

where  n = reflection  (visual  albedo) 

<p  - sphere  illuminance 

<x  - illumination  phase  angle  (apex  at  sphere) 

Clearly,  the  mean  brightness  of  the  sphere  due  to  direct  solar  illuminance 
is  simply  through  interpretation  of  cf>  and  a as  solar  values. 

In  addition,  however,  we  do  not  ignore  the  illumination  of  the 
sphere  by  sky  and  ground,  and  an  account  of  these  contributions  is  also  included. 
Only  for  the  purpose  of  this  calculation  do  we  assume  the  sky  and  ground  to  be 


2.  Russell,  H.N.,  "On  the  Albedo  of  the  Planets  and  Their  Satellites," 
Astrophys.  J.,  43,  3,  pp.  179-196,  April  1916. 
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uniform  hemispherical  illuminators  with  brightness  and  For  horizon- 

'<'7 

tal  viewing,  conservation  and  symmetry  arguments  yield  target  brightness  values 
due  to  the  sky  {B  T 4ACy)  and  ground  &r<i(lD\ 


3 r skv  * z 3 sky 


— B 


( horizontal ) 


For  slant  viewing,  a dependence  on  target  elevation  (€)  angle  was  found  by 
numerical  integration  of  Equation  (1) , and  a good  fit  is  given  by: 


3r’3rMo*iz  0 Sine)  (slant) 


( + GRD  \ 
\ - SKV  ) 


where  the  horizontal  values  {3  T ^0RIZ)  come  directly  from  Equation  (2). 

The  final  step  to  the  determination  of  inherent  target  brightness  is 
to  connect  z?rs/fy,  and  Qr  to  <t>e,  our  selected  common  base.  First,  for  the 
sky,  the  Handbook  of  Geophysics3  conveniently  gives  ground  illumination 
data  broken  down  into  direct  solarCi^)  and  sky  only  (/s<y ) components  as  a 
function  of  sun  elevation.  We  express  these  in  a theoretical  form  as: 


^ o ~~  ^ o ® 1 H h o 


^ 3 Sky 


and  therefore. 
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Using  the  data,  numerical  evaluation  of  Equation  (5)  yields  that,  for  all 
but  the  very  smallest  sun  elevations,  this  ratio  is  almost  constant  with  a 


value  of  0.045.  Thus: 


&SK*  ' 0.  0*5  <p0 


3.  United  States  Air  Force,  Handbook  of  Geophysics,  MacMillan,  1961. 


The  ground  is  now  taken  to  be  a diffuse  reflector  with  albedo  . Both  direct 
sunlight  and  the  sky  contribute: 
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The  contribution  to  inherent  target  brightness  by  the  diffuse  sphere  model 
is  now  complete. 


Section  3 
FLIGHT  TEST 


On  25  February  1976,  a flight  test  was  conducted  in  a desert  region 
by  MICOM/USAADS . The  purpose  of  the  test  was,  in  part,  to  collect  the 
data  necessary  for  registration  of  the  diffuse  sphere  model. 


3.1 


TEST  DESCRIPTION 


The  test  was  constructed  as  four  missions,  each  of  which  involved 
several  passes  of  a low  altitude  jet  fighter  over  a data  collection  site. 

The  passes  covered  a variety  of  approach  azimuths  in  order  to  provide  variable 
sun  angle.  Although  the  passes  were  close  to  radial,  the  natural  offsets 
were  sufficiently  large  that  useful  data  were  obtained  for  aspects  other 
than  head-on.  Overall,  the  four  missions  contained  45  passes. 

The  resulting  test  data  useful  to  this  study  were  of  three  types: 

(1)  Photographic  imagery  of  aircraft  and  luminance  calibration  panels 

(2)  Radar  coordinates  of  aircraft 

(3)  Camera  operator's  log 

The  prime  data,  the  photographic  imagery,  was  collected  on  Plus  X film  with 
a Calspan-owned  system  consisting  of  a Cine-Kodak  Special  II  16  mm  camera 
with  a Century  Tele-Athenar  500  mm  lens.  A Tiffen  #11  filter  was  used  to 
secure  a photopic  sensitivity.  The  frame  rate  (for  this  test)  was  set  at 
one  per  second,  and  tracking  was  facilitated  by  a Redfield  3X-9X  zoom 
riflescope  boresighted  with  the  camera.  On  each  pass,  photographic  recording 
was  initiated  as  soon  as  the  aircraft  was  acquired  by  the  camera  crew. 

Usually,  the  recording  ended  near  crossover,  but  occasionally,  outbound 
frames  were  collected. 

In  addition  to  the  photography,  radar  coordinates  vs  time-of-day 
were  provided  during  each  pass  at  a rate  of  one  per  second.  Although  the 
camera  and  radar  data  rates  were  the  same,  they  were  not  actually  synchronized. 
The  precision  of  the  radar  data  was  more  than  sufficient  for  the  needs  of  the 
present  study. 
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The  camera  log  contained  the  detailed  flow  of  the  test.  Included 
were  camera  calibration  data,  definitions  of  missions  and  passes,  and  comments 
on  the  passes  wnere  appropiate  such  as  camera  start  and  step  times  and 
estimated  number  of  frames  collected. 

3.2  PHOTOMETRIC  CALIBRATION 

Special  emphasis  was  placed  on  calibration  of  the  camera  system  to 
assure  that  absolute  luminance  values  could  be  obtained  from  the  film.  Prior 
to  the  test,  a series  of  exposures  was  taken  over  the  gamut  of  exposure  times 
and  stops  available  on  the  camera.  In  support  of  these  exposures,  source 
luminance  measurements  were  taken  with  a Pritchard  model  1970  photometer. 

The  resulting  film  densities  could  now  be  related  to  source  luminance, 
exposure  time,  and  f/number.  This  relationship,  together  with  a criterion 
that  sky  brightness  yields  film  densities  in  the  range  of  1.0-1. 5,  became  a 
rule  of  thumb  for  exposure  time  and  f/number  setting.  Accordingly,  at  the 
time  of  the  test,  sky  brightness  measurements  near  the  horizon  indicated 
1/70-second  exposure  at  f/22.  These  settings  were  used  throughout  the  test. 

The  calibration  discussed  above  was  used  only  for  camera  setup. 

The  absolute  connection  between  film  density  and  source  luminance  was  obtained 
by  photographing  two  "field  step  wedges"  located  a short  distance  away  from 
the  camera.  One  wedge  faced  south  and  the  other  north  to  obtain  increased 
brightness  dynamic  range  according  to  sun  angle.  Each  wedge  consisted  of 
eight  sheets  of  matte  paper  of  differing  shades  of  gray.  The  sheets,  arranged 
in  a 2 x 4 array,  had  approximate  reflectivities  of  1.00,  0.70,  0.48,  0.35, 

0.24,  0.13,  0.07,  and  0.03.  Both  wedges  were  photographed  at  times  of 
opportunity  (between  missions)  in  conjunction  with  remote  luminance  measure- 
ments with  the  Pritchard  photometer  of  each  of  the  eight  steps.  These  measure- 
ments directly  provided  absolute  film  calibration,  and  it  is  noted  that  knowledge 
of  reflectivity,  such  as  the  values  quoted  above,  is  unnecessary. 
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Section  4 
DATA  ANALYSIS 


4.1  INTRODUCTION 

With  the  exception  of  the  details  of  this  diffuse  sphere  model 
(Section  2) , all  data  reduction  and  analytic  procedures  are  described  in  this 
section.  This  includes  reduction  of  the  photographic  data  (together  with 
merger  of  appropriate  extracts  from  the  radar  data)  and  registration  of  the 
diffuse  sphere  model. 

Inherent  in  the  model  registration  is  an  account  of  atmospheric 
effects.  The  diffuse  sphere  model  treats  only  the  intrinsic  brightness  of 
the  aircraft,  whereas  the  photographic  data  represent  apparent  brightness 
according  to  camera  position.  Along  the  path  between  the  two,  both  atmospheric 
attenuation  and  inscattering  of  light  occur,  and,  to  make  the  connection,  a 
description  of  these  effects  is  included.  Since  the  purpose  of  the  investi- 
gation is  a test  of  the  diffuse  sphere  model,  as  much  detail  as  possible  was 
included  on  atmospherics  to  eliminate  it  as  an  issue.  Accordingly,  some  of  the 
following  subsections  become  rather  involved  with  atmospherics,  but  the  reader 
should  not  conclude  from  this  that  atmospheric  effects  are  correspondingly 
important.  Actually,  visibility  during  the  test  was  excellent,  and  the  photo- 
graphic ranges  were  short.  Although  atmospherics  did  affect  luminance  values, 
it  was  not  a major  factor. 

4.2  ANGULAR  DATA 

Both  the  diffuse  sphere  and  the  atmospheric  models  require  certain 
angular  data  as  inputs.  As  a matter  of  convenience,  for  each  data  point  in  the 
radar  record,  aircraft  azimuth  and  elevation,  sun  azimuth  and  elevation,  and 
solar  phase  angle  (angle  between  the  sun  and  camera  with  apex  at  the  aircraft) 
were  computed  and  made  available  as  a data-set.  Later  in  the  analysis,  extracts 
from  these  data  are  merged  with  the  photographic  data. 
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4.3 


DENSITOMTRY 


The  first  step  in  the  analysis  of  the  film  data  is  the  reduction  of 
the  photographic  images  of  aircraft  and  sky  to  suitable  numerical  form.  To 
accomplish  this,  use  was  made  of  a two-dimensional  scanning  densitometer 
manufactured  by  Spatial  Data  Systems.  This  instrument  consists  of  a closed- 
circuit  television  system  with  the  ability  to  provide  values  proportional 

■ 

to  total  areas  within  the  scanned  field  that  lie  within  predetermined  brightness 
values.  Thus,  by  providing  the  TV  sensor  with  a light  table  and  appropriate 
optics,  individual  16  mm  frames  were  scanned  and  the  image  dichotomized  into 

i 

a sequence  of  numbers  representing  areas  within  the  frame  corresponding  to 
neighboring  density  increments  of  the  film. 

The  brightest  level,  the  darkest  level,  and  the  number  of  steps 
* between  are  degrees-of-freedom  of  the  densitometer,  and  their  selection 

entailed  examination  of  the  film  data  in  advance  of  the  final  measurements. 

It  was  therefore  determined  that  one  of  the  frames  containing  a brighter 
(south-facing)  image  of  the  eight-step  field  wedge  was  satisfactory  for 
calibration,  and  a procedure  for  "locking  on"  the  wedge  with  the  brightest 
and  darkest  levels  was  established. 

By  trial,  it  was  found  that  16  bands  or  density  increments  between 
the  two  limits  could  be  included  for  measurement  without  undue  hardship.  In 
' regard  to  the  increments,  two  comments  are  appropriate.  First,  the  densitometer 

allows  the  selection  of  the  scaling  of  these  increments  to  be  linear  (approxim- 
ately) with  either  film  density,  which  varies  as  the  log  of  film  transmission, 
or  transmission  itself.  Since  our  interest  will  be  in  averaging  over  scene 
luminance,  which  relates  closely  to  film  transmission,  the  latter  mode  was 
employed.  Secondly,  a valuable  adjunct  to  the  densitometer  is  a monitor  which 
displays  the  scanned  image,  with  the  sequence  of  transmission  increments 
appearing  as  individual  false  colors.  This  is  especially  useful  for  identi- 
fication of  areas  within  the  aircraft  silhouette  that  fall  within  the  same 
brightness  increment  as  the  sky  (or  terrain)  present  in  the  frame.  In  these 
cases,  since  the  regions  involved  external  to  the  aircraft  image  consume  large 
areas  of  the  frame,  the  area  measures  for  those  increments  are  contaminated 
beyond  use.  In  all  cases  where  this  occurred,  visual  estimates  were  made  of  the 
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percent  of  silhouette  area  for  each  increment  involved.  The  incorporation  of 
these  estimates  into  the  analysis  will  be  discussed  later. 

} 

The  calibration  was  completed  by  meshing  the  photometer-measured 
luminance  values  for  the  eight-step  wedge  and  the  16  increments  to  be  used  for 
data  reduction.*  The  results  of  this  yielded  a luminance  schedule  for  the 
increments  as  follows: 


Increment 

Luminance  (ft-L) 

Increment 

Luminance  (ft-L) 

1 

50 

9 

1600 

2 

180 

10 

1890 

3 

270 

11 

2190 

4 

370 

12 

2550 

5 

500 

13 

2980 

6 

700 

14 

3730 

7 

990 

15 

4500 

8 

1280 

16 

5260 

The  luminance  values  represent  the  centers  of  the  increments. 

For  each  pass  containing  usable  imagery,  one  of  the  frames  near 
crossover  was  selected  as  number  one  in  an  ad-hoc  frame-numbering  sequence,  with 
frame  number  increasing  with  successive  earl ier  frames,  i.e. , backward  on  the 
film  strip.  Since  the  imagery  was  usually  of  the  aircraft  while  inbound, 
increasing  numbers  meant  greater  range  and  correspondingly  smaller  image  size. 
Successive  frames  were  analyzed  in  increasing  number  until  either  camera 
acquisition  was  reached,  or  resolution  became  a problem.  The  results  of  this 
analysis  were  punched  on  IBM  cards  to  be  available  as  data  for  the  next  data 
reduction  stage.  Each  frame  analyzed  was  represented  by  one  card  containing 
pass  number,  frame  number,  area  measures  (or  % areas)  for  the  16  increments, 
increment  number  for  the  sky  near  the  aircraft,  and  the  increment  number  for 
terrain  if  present  in  the  frame.  A printout  of  these  data  is  given  in  Appendix  A. 

★ 

For  calibration  only,  32  increments  were  actually  used  to  obtain  increased 
luminance  resolution. 


w 
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4.4  MEAN  AIRCRAFT  LUMINANCE 

The  diffuse  sphere  model  described  in  Section  2 predicts  average 
luminance  of  the  spherical  silhouette.  Since  we  have  no  specific  interest 
in  the  luminance  structure,  the  mean  of  the  data  can  be  computed  now.  Although 
the  data  represent  aircraft  luminance  apparent  to  the  film  and  the  model 
represents  intrinsic  luminance  of  the  aircraft,  corrections  for  atmospherics 
and  camera  flare  (anticipated  here)  can  be  treated  independently.  That  is, 
the  triad  of  operations: 

(1)  computation  of  silhouette  mean  luminance 

(2)  correction  for  atmospherics 

(3)  correction  for  flare 

are  all  processes  obeying  the  principle  of  linear  superposition  and  can  be 
freely  transposed. 

An  expression  for  the  mean,  accounting  for  estimated  percentages  in  the 
absence  of  measured  areas,  is  now  derived.  The  mean  luminance  (5)  is  written: 

21  aiBi  fs'i 


where : 


A » 21  ai 

4-  1 


The  summations  are,  of  course,  over  the  N non-zero  increments  of  the  data  in 
Appendix  A.  For  our  purpose  here,  let  the  first  M components  of  the  summations 
be  the  measured  areas,  and  the  components  M+l  to  N,  estimated  percentages. 
Equation  (8)  is  now  expressed: 


3 = -j  (ll  «4  Bi  + aJ8J 

^ \t'1  y 
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Eliminating  A between  Equations  (11)  and  (12)  yields: 


1 ~ Z 

j-M  + 1 , 


z * 

-i'-c  1 

Af 

zx 

«*/ 


* £ 


which  is  the  final  expression  to  be  used. 

Before  leaving  the  concept  of  incorporating  estimated  percentages, 
a comment  is  in  order  regarding  the  stability  of  B on  errors  in  estimation. 

We  note  that  for  a small  unknown  a^ , the  estimate  f ^ has  high  relative  error. 
However,  since  the  contribution  to  B is  correspondingly  small,  no  instability 
occurs.  On  the  other  hand,  as  a^  and  f.  become  large,  the  relative  error  of 
f j becomes  small.  Although  absolute  errors  in  f^  (when  a^  is  large)  cause 
large  relative  errors  in  the  contributions  of  the  measured  a^'s,  they  are 
themselves  small,  and  again,  no  instability  occurs.  On  one  frame  (pass  27, 
frame  12),  it  was  necessary  to  estimate  all  non-zero  increments,  but  the 
resulting  mean  luminance  was  nonetheless  acceptable. 


At  this  point,  some  emending  of  the  film  data  was  necessary.  Due 
to  low  offset  values,  high  angular  rates  were  possible  near  crossover,  thereby 
indicating  a high  time  resolution  requirement  for  camera/radar  data  synchro- 
nization. Since  it  was  indicated  in  the  camera  log  that  camera  times  near  the 
end  of  the  passes  were  of  questionable  accuracy,  they  were  not  employed. 
Instead,  the  film  was  previewed  to  locate  crossover,  a matter  of  visual  aspect 
of  the  aircraft  engine  exhaust  port.  Since  the  frame  nearest  crossover  became 
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the  key  to  synchronization  with  the  radar  data,  all  passes  that  did  not  contain 
crossover  imagery  were  deleted  from  further  analysis.  In  addition,  anomolous 
low  luminance  values  were  noted  from  the  densitometry  of  passes  36  through  42, 
and  these  were  deleted  as  well. 

The  surviving  data  were  analyzed  according  to  Equation  (13),  and, 
as  a matter  of  general  interest,  apparent  contrast  (C)  of  the  aircraft  according 


B - Bc 


was  computed  as  well.  Bs  is  selected  from  the  schedule  in  section  4.3  according 
to  the  sky  brightness  increment  number.  By  means  of  identifying  the  crossover 
frame  within  the  frame  scale  for  each  pass,  the  radar  and  camera  data  were 
synchronized.  This  allowed  B and  the  associated  geometry  parameters  necessary 
for  model  registration  to  be  merged.  The  merged  data  were  punched  on  IBM  cards, 
with  one  card  again  representing  one  frame,  as  follows: 

(1)  Pass  number 

(2)  Frame  number 

(3)  Number  identifying  radar  data  point 

(4)  Slant  range  (meters) 

(5)  Aircraft  elevation  angle  (degrees) 

(6)  Solar  elevation  angle  (degrees) 

(7)  Solar  phase  angle  (degrees) 

(8)  Sky  luminance  (ft-L) 

(9)  B,  mean  apparent  aircraft  luminance  (ft-L) 

A printed  list  of  this  data  set  is  included  in  Appendix  B.  The  printed  version 
includes  the  resulting  apparent  contrast  values  as  well. 

In  anticipation  that  solar  phase  angle  is  the  dominant  independent 
variable  of  the  diffuse  sphere  model,  the  mean  illuminance  data  were  simply 
plotted  against  phase  angle.  The  result  is  shown  in  Figure  1,  and  three 
observations  are  made: 
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luick  Look  at  Mean  Luminance  Values 


(1)  At  angles  less  than  90°,  the  decreasing  luminance 
with  increasing  angle  is  consistent  with  the  model, 
although  the  slope  is  less  than  anticipated. 

(2)  At  angles  greater  than  90°,  the  increasing  luminance 
with  increasing  angle  is  in  direct  opposition  to  the 
model  since  the  model  luminance  is  monotonic  in  phase 
angle. 

(3)  In  the  mid-range  of  phase  angle,  three  data  points  lie 
significantly  above  the  main  trend. 

In  regard  to  (3),  it  was  suspected  that  solar  glints  were  responsible.  Glint 
is  loosely  defined  here  as  those  singularly  intense  specular  reflections  from 
surfaces  that  are  close  to  being  either  flat  or  conical.  (It  is  noted  that 
the  skin  of  the  aircraft  is  unpainted  metal  and  specular  reflections  can  be 
expected.)  The  three  frames  involved  (pass  3,  frame  2 and  pass  23,  frames  4 
and  5)  did  indeed  contain  significant  glint  and  were  not  included  in  the  model 
registration.  However,  the  occurrence  of  these  glints  will  not  be  ignored  in 
the  final  evaluation  of  the  model. 

The  departure  indicated  in  (2)  has  major  significance.  If  it  does 
not  represent  a failure  on  the  part  of  the  diffuse  sphere  model,  it  must  be 
due  to  a mechanism  that  causes  aircraft  luminance  apparent  at  the  film  to 
differ  from  inherent  aircraft  luminance  in  a fashion  consistent  with  the  data, 
and  two  candidates  can  be  identified: 

(1)  Camera  flare  from  direct  solar  irradiance 

(2)  Strong  forward  scattering  by  the  atmosphere 
(i.e.,  Mie  scattering  by  haze) 

Thus,  it  is  critical  to  the  resolution  of  the  question  of  the  departure  at 
large  phase  angles  that  both  these  mechanisms  be  accounted  for  at  the  time  of 
model  registration,  and  this  forms  the  subject  of  the  following  section. 
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4.5  APPARENT/ INTRINSIC  LUMINANCE  CONNECTION 

We  now  address  the  connection  between  apparent  and  intrinsic  aircraft 
luminances  in  terms  of  camera  flare  and  atmospherics  that  have  been  identified 
as  accountable  mechanisms  prior  to  model  registration. 

For  flare,  the  camera  system  underwent  laboratory  measurements  to 
determine  film  plane  irradiance  as  a function  of  phase  angle.  Scaling  to  the 
flight  test  illumination  levels  was  accomplished  through  the  luminance  value 
measured  for  one  of  the  steps  in  the  field  calibration  wedge  and,  accounting 
for  sun  angle,  an  equivalent  laboratory  simulation.  The  measurements  show 
flare  can  be  neglected  for  angles  less  than  165°.  For  angles  increasing  above 
165°,  flare  increases  very  rapidly,  reaching  580  ft-L  at  172°.  The  accuracy 
of  the  laboratory  data,  both  in  angle  and  source  collimation,  precluded  attempts 
to  provide  a phase-angle-dependent  flare  correction.  Accordingly,  deletion  of 
data  at  these  high  angles  (passes  44  and  45)  from  model  registration  eliminated 
camera  flare  as  a factor. 

The  remainder  of  the  discussion  here  presents  a mathematical  system 
for  handling  the  optical  effects  of  the  atmosphere.  Although  the  system  may 
appear  relatively  involved,  it  does  not  necessarily  mean  that  atmospherics  are 
important.  It  only  means  the  descriptive  tools  were  readily  available,  and 
to  identify  and  sift  out  those  that  would  be  unimportant  seemed  pointless.  As 
a matter  of  hindsight,  atmospherics  did  not  turn  out  to  be  the  prime  factor 
in  model  registration. 

The  atmosphere  has  two  distinct  optical  effects  beyond  the  illumina- 
tion scheme  discussed  in  Section  2.  First,  propagation  of  luminance  follows  a 
transport  equation  involving  atmospheric  parameters  for  points  along  the  path 
between  the  aircraft  and  camera,  and,  second,  the  low-altitude  illuminance  <#> 
[which  appears  in  Equation  (1)  as  well  as  below]  is  reduced  according  to  solar 
elevation  angle  through  the  total  air  mass  above  the  earth.* 


The  reader  not  specifically  interested  in  deta;ls  of  atmospheric  effects  can 
proceed  to  the  last  sentence  before  section  4.6,  and  continue. 
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The  propagation  will  be  handled  here  by  a standard  form  of  the 
transport  equation  known  as  Koschmieder ' s Equation:* 


Bo  G 


-T 


3oo(>  - e~T) 


(15) 


where: 


~ o 

B 

r 

B oo 


- inherent  luminance  of  aircraft 
= aircraft  luminance  apparent  at  the  camera 
= optical  thickness  of  aircraft/camera  path 
= atmospheric  "equilibrium"  luminance 


The  conditions  of  validity  are  that,  over  the  path  in  question,  both  the  solar 
irradiance  and  the  shape  of  the  scattering  coefficient  (phase  function)  are 
stationary.  Since  the  aircraft  elevation  is  low  (usually  a few  hundred  meters) 
and  atmospherics  is  not  expected  to  be  a major  factor,  variability  of  atmos- 
pheric properties  over  the  path  will  not  be  addressed. 

In  Equation  (15) , the  optical  thickness  is  simply: 


r = jSR 


where:  /3  = attenuation  coefficient  (m  *) 

R = aircraft/camera  slant  range  (m) 

with /3  in  turn  given  by  the  definition  of  meteorological  range  (V) 


(16) 


A = 


brx  5~0 


(17) 


The  terminology  used  here  and  in  the  equations  that  follow  are  standard. 

For  further  details,  the  equation  is  derived  from  a more  general  form  in 
Appendix  C,  which  includes  definitions  and  validity  conditions.  Beyond  this, 
any  number  of  texts  can  be  found.  Middleton  is  especially  lucid. 
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The  equilibrium  luminance,  , is  given  by  the  spherical  integral: 

/^oo  = $ ^ &su*.  d .n. 


A 

where:  k = phase  function  at  path  point 

&sux  = spherical  luminance  function  surrounding  path  point 

and  the  integration  is  over  4)T  steradians  about  the  point.  To  obtain  the 
phase  function,  we  recognize  two  independent  physical  scattering  mechanisms: 
scattering  by  atmospheric  density  fluctuations  on  a molecular  scale  (Rayleigh 
scattering)  and  scattering  by  water  haze  (Mie  scattering) . The  form  for 

A 4 

Rayleigh  scattering  (kR)  is  well-known  : 

**  = jh  * cos 2T)  (19) 

where  y is  the  scattering  angle.  (For  the  direct  sunlight  component,  y is 
simply  T - a,  the  phase  angle.)  The  phase  function  is  very  much  different  for 
Mie  scattering  (k^i . Data  given  by  Deirmendj ian ^ can  be  expressed  by  a 
simp! e model : 

{,.  100-73  - ^r-oW  (20) 

A 

with  y expressed  in  radians.  A plot  of  kM  vs  y in  degrees  is  given  in 
Figure  2 , and  very  strong  scattering  in  the  forward  direction  is  evident. 

This  is  precisely  the  effect  described  earlier  as  being  important  that  it 
not  be  overlooked  as  a possible  explanation  for  the  trend  of  the  data  at  high 
solar  phase  angles. 


4.  S.  Chandrasekhar,  Radiative  Transfer,  Dover  Press  (1960) 

5.  D.  Deirmendjian,  Electromagnetic  Scattering  on  Spherical  Polydispersions, 


Elsevier  Publishing  Co.  (1969) 
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With  the  individual  scattering  mechanism  described,  the  resultant 

A 

phase  function  k required  by  Equation  (18)  is  given  by  the  weighted  super- 
position of  the  components: 

k - (21) 

/3 


where  /5R  and  /5„  are  the  individual  Rayleigh  and  Mie  scattering  coefficients 


which  satisfy: 


A - Ax  * A*  (22) 

The  Rayleigh  scattering  coefficient,^^,  is  a constant,  and  we  will  give  its 
value  below  along  with  other  constants.  The  Mie  coefficient,  which  is 
proportional  to  haze  concentration,  is  a variable,  and  its  value  is  found 
from  Equation  (22)  in  conjunction  with  Equation  (17) . 

We  now  turn  our  attention  to  the  completion  of  the  spherical  inte- 
gration of  Equation  (18) . Through  a piecewise  breakdown  of  the  integrand  and 
a few  approximations,  a simple  algebraic  expression  for  can  be  obtained. 

First,  the  surround  luminance,  BCIID,  can  be  written  as  the  sum  of  direct 

bUK 

sunlight  B0  and  a residual  diffuse  part  Bpjp: 


BSUR  B©  + bdif 


and  Equation  (18)  becomes: 


= j>  k 3Q  dn.  + <f>  k &Dn.  dj-i. 


Since  B0  is  bounded  to  the  solar  disc,  which  has  a subtense  only  1/2  in 

A A 

diameter,  k is  constant  and  assumes  a value  (kQ)  according  to  the  solar  phase 
angle.  We  can  now  write: 


j)kBQd^-  - k0  B0  ddi  = kQ  0e 
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where  is  the  direct  solar  irradiance.  In  regard  to  the  diffuse  part  of 
Equation  (24),  a further  expansion  is  obtained  through  the  separate  scattering 

, A 

mechanisms.  Elimination  of  k between  Equation  (21)  and  the  diffuse  term 
yields : 

$ k Bp, p cIj~l  - -g-  (^/3„  Bdif  d-H.  +/?%  (f>  Bj>if  ^ (26) 

A A 

Now,  since  ^(7-)  sin  y is  concentrated  at  small  f as  well  as  kM(y) , the  main 
contribution  to  the  first  integral  in  Equation  (26)  is  from  BDIp  in  the 
vicinity  of  y = 0.  We  therefore  make  the  approximation  that  is  a 

constant  equal  to  the  sky  brightness,  , extracted  from  the  film.  Thus: 

$ - &s  (27) 


The  final  integral,  the  diffuse  term  of  Equation  (26),  possesses  an 
integrand  whose  dependence  on  y is  very  weak.  Ignoring  the  ground  as  a con- 
tributor, Bpjp  is  replaced  with  a representative  mean  value  for  the  sky  and, 
for  this  purpose,  the  value  Bg^y  [Equation  (6),  Section  2]  is  used.  Therefore, 
we  have: 

$ &PIF  = ~2  &SKY  (28) 

and  the  integration  to  yield  BMis  complete.  Collecting  all  the  pieces. 
Equation  (18)  is  finally  written: 


3 


00 


0.0225 


Am 


(29) 


What  now  remains,  to  complete  the  connection  between  apparent  and 
intrinsic  aircraft  luminance,  is  a procedure  for  generating^,  and  the 
specification  of  some  left-over  constants.  The  parameter  used  both  here 
and  in  Section  2,  is  direct  solar  illuminance  (perpendicular  to  sun  direction) 
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at  ground  level.  As  the  elevation  angle  of  the  sun  is  reduced,  (p&  decreases 
due  to  an  increased  path  through  the  atmosphere  according  to: 

- r* 

4"  ^oe  (30) 

where  0QO  = solar  illuminance  outside  atmosphere 

T* = vertical  optical  thickness  of  air  mass 
hQ  = solar  elevation  angle 

6 

From  an  atmospheric  model  by  Elterman  , T can  be  found: 


= °-f 47  + fa  (31) 

The  parameter is  the  Mie  scattering  coefficient  at  an  altitude  of  5000  meters, 
a pivot  point  in  this  model.  The  value  for  /Ss  , as  well  as  the  value  for /3K 
needed  in  earlier  equations,  are  given  by  Elterman  as: 

fis  = 5.02  x 10  ^ meters" 1 

-5  -1  (32> 

ft*  = 1.162  x 10  meters 

The  connection  between  apparent  and  intrinsic  aircraft  luminance 
is  now  complete  to  within  the  two  parameters,  V (meteorological  range)  and  0OO 
(a  number  for  solar  illuminance). 


4.6  MODEL  REGISTRATION 

The  atmosphere  description,  together  with  the  diffuse  sphere  model 
of  Section  2,  form  the  total  model  for  registration  on  the  data.  Overall, 
four  parameters  remain  unspecified: 


6.  L.  Elterman,  "UV,  Visible,  and  IR  Attenuation  for  Altitudes  to  50  KM,  1968", 
AF  Cambridge  Research  Laboratories,  AFCRL-68-0153  (April  1968). 
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(1)  a , aircraft  reflectance  (albedo) 

(2)  /i  , ground  reflectance  (albedo) 

(3)  d>00,  solar  illuminance  external  to  atmosphere 

(4)  V , meterological  range  (meters) 

1 

Reflectance  values  for  a variety  of  terrains  are  given  by  Gorden  and  Church  . 
From  these,  a value  of  0.22  ("Dirt  - flat  desert  road  freshly  bulldozed")  was 
assigned  to  . The  remaining  three  parameters  were  varied  in  initial 
attempts  to  effect  a model/data  fit. 

The  criterion  for  registration  was  to  minimize  the  total  square 
error  in  apparent  aircraft  luminance  by  means  of  a global  search  in  the  space 
of  the  adjustable  parameters.  Such  a search  over  tne  space  of  , 0Oo,  and  V 
did  not  converge  to  a reasonable  solution.  Apparently,  the  Mie  scattering 
mechanism,  although  having  a proper  phase  angle  dependence  as  noted  earlier, 
was  too  weak  at  normal  meteorological  ranges;  and,  to  compensate,  the 
registration  increased , thereby  driving  V to  an  extremely  small  value. 

For  proper  registration,  it  became  necessary  to  specify  both  V and yt 
and  to  leave  <pOQ  as  the  only  variable.  A value  of  40  km  was  selected  to  be  as 
low  as  the  writer  was  willing  to  use  in  light  of  the  observation  that  visibility 
during  the  test  was  "very  high."  Should  atmospherics  remain  an  issue,  further 
consideration  would  be  given  to  V as  required. 

The  surface  of  the  aircraft  was  primarily  (but  not  entirely)  unpadnted 
metal,  and  a value  of  0.7  was  chosen  for  /i  • This  is  moderately  less  than  the 
visual  albedo  for  aluminum.  The  importance  of  the  reflectance  value,  again,  is 
related  to  the  importance  of  the  atmospherics.  An  error  in  reflectance 
will  in  part  be  compensated  for  by  <POQ  during  registration  since  intrinsic 
aircraft  luminance  is  proportional  to  the  product  AtP^c-  The  compensation  is 
not  complete  because  affects  apparent  luminance  through  inscattering, 
whereas yi  does  not.  Thus,  if  inscattered  light  between  the  aircraft  and  camera 
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7. 


J.I.  Gorden  and  P.V.  Church,  "Directional  Luminous  Reflectances  of  Objects 
and  Backgrounds  for  a Moderately  High  Sun,"  Contained  in  S.Q.  Duntley  et  al . , 
"Visibility,"  University  of  California  Scripps  Institute  of  Oceanography. 
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is  not  important,  neither  is  aircraft  reflectance.  We  shall  treat /t  in  the 
same  fashion  as  V,  in  that  further  consideration  of  its  value  will  be  given 
if  necessary. 


Model  registration  was  completed  through  adjustment  of  the  lone 
parameter  00o.  The  completed  model  was  then  exercised  using  as  input  all 
necessary  data  corresponding  to  each  measured  apparent  luminance  value,  to 
yield  predicted  apparent  luminances  on  a one-for-one  basis.  The  resulting 
data/model  luminance  pairs  are  plotted  in  (a)  and  (b)  of  Figure  3 against 
solar  phase  angle.  (The  division  into  two  plots  was  only  to  facilitate 
presentation  of  the  large  number  of  points;  no  distinction  between  the  missions 
is  intended.)  From  these  plots,  the  following  observations  are  made: 

(1)  In  the  range  of  0-90°,  the  trend  of  decreasing  luminance 
with  phase  angle  is  present  in  the  data  but  is  weaker 
than  the  model  predictions. 

(2)  At  angles  greater  than  90°,  significant  departures 
occur  between  the  data  and  the  model.  Since  forward 
scattering  by  haze  was  not  significant,  the  departures 
are  considered  real,  in  that  they  are  present  in  the 
intrinsic  luminance  as  well  as  the  apparent  luminance. 


(3)  The  overall  departure  of  the  model  from  the  data  is  much 
greater  than  deviations  within  the  model  at  any  given 
phase  angle.  Since  these  deviations  contain  atmospheric 
effects  along  the  slant  path  between  the  aircraft  and 
camera,  atmospherics  can  be  regarded  as  unimportant.  Thus, 
the  selection  of  values  for  meteorological  range  and 
aircraft  reflectance  need  not  be  redressed. 


(4)  In  regard  to  deviations  within  the  model,  two  model  points 

lie  below  2000  ft-L  at  small  phase  angles  on  (b)  of  Figure  3. 
These  two  points  occurred  on  pass  35  with  a solar  elevation 
angle  of  only  19°,  and  significant  downward  correction  in 
aircraft  illuminance  resulted  from  attenuation  by  the  total 
air  mass. 
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4.7  THE  GLINT  HYPOTHESIS 

This  section  concludes  with  the  suggestion  that  glint  is  likely 
responsible  for  the  principal  model/data  departure.  Although  inconclusive, 
the  question  of  whether  or  not  the  geometry  is  correct  for  glints  to  occur 
can  be  answered  approximately.  To  do  this,  a sufficient  first  order  descrip- 
tion of  the  aircraft  is  as  follows: 

(1)  Wings  - horizontal  flats 

(2)  Tail  - vertical  and  horizontal  flats 

(3)  Fuselage  - cylinder 

Although  the  wings  are  an  unlikely  source  of  glint,  all  surfaces  as  represented 
above  have  surface  normals  at  90°  to  the  forward  axis  of  the  aircraft.  Since 
specular  reflections  occur  when  the  solar  phase  angle  is  bisected  by  the 
surface  normal,  the  data  we  have  can  be  put  to  a test.  Thus,  for  all  the 
data,  the  angle  between  the  phase  angle  bisector  and  the  forward  direction  of 
the  aircraft  was  computed  and  plotted  against  phase  angle  itself  (Figure  4) . 

The  encircled  points  are  those  in  Figure  3 having  the  large  departures  from  model 
prediction.  They  primarily  fall  in  the  range  of  100-110°  aft  of  the  forward 
direction,  with  some  at  120°.  Thus,  there  is  an  excellent  agreement  between 
the  phase  angle  bisector  and  normals  to  aircraft  surfaces  that  we  would  expect 
to  be  responsible  for  glint.  This  is  regarded  as  strong  evidence  that  the 
primary  departure  between  the  data  and  the  model  is  due  to  glint. 
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Figure  4 Connection  Between  Phase  Angle  Bisector  Direction  and  Phase  Angle 
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Section  5 
CONCLUSIONS 

The  prime  result  of  this  study  was  the  registration  of  a diffusely 
reflecting  sphere  model  for  predicting  aircraft  luminance  on  field  test 
data,  and  the  resulting  comparison  was  shown  in  Figure  3.  The  question  of 
whether  or  not  the  model  is  valid  or  appropriate  requires  a highly  qualified 
answer.  It  depends  on  the  accuracy  required  in  application,  the  nature  of 
the  use  of  the  results,  and  the  availability  or  appropriateness  of  alter- 
nate models. 

For  example,  the  model  should  probably  not  be  used  when  sun  angle 
is  a critical  issue.  On  the  other  hand,  applications  that  only  require  a 
measure  of  variability  due  to  sun  angle  may  find  the  model  useful.  The 
answer  as  to  whether  or  not  the  model  is  appropriate  must  come  from  the  user. 

Probably  the  greatest  limitation  of  the  model  is  the  availability 
of  significantly  better  models  not  much  more  complex.  The  crude  geometrical 
representation  of  the  aircraft  given  in  Section  4.7,  or  something  akin  to  it, 
could  be  the  basis  of  a different  model  resulting  in  a much  better  fit  to 
the  field  data.  The  fact  that  the  aircraft  used  in  the  test  was  neither 
diffuse  nor  a sphere  is  not  a fault  of  the  model,  since  its  use  was  strictly 
empirical.  However,  by  including  a little  attention  to  the  actual  surface 
shape,  and  the  reflective  properties,  improvements  probably  will  occur.  It 
is  noted  that  optical  reflection  signatures  for  diffuse  and  specular  flats, 
cones,  and  cylinders,  and  compound  surfaces  such  as  spheres  are  known,  and 
a replacement  for  the  diffuse  reflecting  sphere  seems  very  close  at  hand. 
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Appendix  A 
DENSITOMETRY  DATA 


Included  here  are  the  raw  data  for  aircraft  brightness  resulting 
from  the  densitometery  analysis  of  the  flight  test  imagery.  The  first  two 
columns  contain  the  pass  number  and  camera  frame  number,  respectively  (P/F). 

The  next  16  columns  contain  the  actual  areas  of  the  aircraft  image  (arbitrary 
units)  that  fall  within  each  of  the  16  brightness  increments  employed  in  the 
dichotomy.  For  those  increments  that  could  not  be  measured  due  to  contamination 
from  regions  external  to  the  aircraft  silhouette,  percentage  estimates  of  total 
silhouette  area  were  made  and  labeled  negative  for  identification.  The  last 
two  columns  are  the  brightness  bin  numbers  for  the  sky  (S)  and  typical  terrain 
(T)  when  present  in  the  frame. 

The  meaning  and  use  of  all  data  in  the  following  table  are  discussed 
in  detail  in  the  main  text. 
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TABLE  A-1.  RAW  DATA  FROM  DENSITOMETER  ANALYSIS 


LUMINANCE  INCREMENT  AREAS  (ARBITRARY  UNITS) 


15 
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TABLE  A-1.  (Cont.)  RAW  DATA  FROM  DENSITOMETER  ANALYSIS 


TABLE  A-1.  (Cont.)  RAW  DATA  FROM  DENSITOMETER  ANALYSIS 


LUMINANCE  INCREMENT  AREAS  (ARBITRARY  UNITS) 

4 5 6 7 8 9 10  11  12  13  14  15 


TABLE  A-1.  (Cont.)  RAW  DATA  FROM  DENSITOMETER  ANALYSIS 


Appendix  B 

DATA  SET  FOR  MODEL  REGISTRATION 


This  appendix  presents  the  merged  data  set  upon  which  model 
registration  was  effected.  In  addition  to  integers  useful  for  identification, 
the  list  includes  all  geometric  and  luminance  values  individually  required 
of  both  the  diffuse  sphere  model  and  the  atmospherics.  The  contrast  values 
were  not  employed  in  the  registration  but  are  included  here  since  they  have 
general  interest. 
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TABLE  B-1.  DATA  ARRAY  FOR  MODEL  REGISTRATION 
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TABLE  B-1.  (Cont.)  DATA  ARRAY  FOR  MODEL  REGISTRATION 


22  6 2363  1367.  3.  33.  **.  2930.  2370.  -0.20 

22  5 236*  1369.  *.  38.  *7.  2930.  231*.  -0.22 

22  3 2366  1053.  5.  33.  5b.  2930.  22*9.  -0.25 

23  5 2*67  1207.  6.  33.  62.  3730.  3027.  -0.19 

23  * 2*69  1015.  7.  33.  77.  2550.  2806.  0.10 

23  3 2*69  353.  9.  33.  69.  2550.  2171.  -0.15 

27  12  2327  2797.  3.  -0.  72.  2550.  1627.  -0.36 

27  11  2829  2559.  3.  *0.  7*.  2550.  1669.  -0.35 

27  10  2929  2323.  3.  *0.  75.  2550.  1536.  -0.*0 

27  9 2330  2039.  3.  *0.  77.  2550.  1593.  -0.38 

27  3 2931  1962.  *.  *0.  9C . 2550.  1555.  -0.39 

27  7 2332  16*0.  *.  *0.  33.  2550.  1500.  -0.*1 

27  6 20  33  1*29.  5.  *0.  08.  2550.  1*69.  -0.*2 

27  5 233*  1235.  6.  *3.  9-.  2190.  1362.  -0.38 

27  * 2835  1066.  7.  *0.  101.  2550.  16*0.  -0.36 

27  2 2337  960.  9.  *0.  12*.  3730.  17u3.  -0.5* 

23  3 2893  2156.  2.  *0.  135.  3730.  2819.  -0.2* 

25  7 289*  1973.  3.  *0.  132.  3730.  2636.  -0.28 

29  6 2395  1819.  3.  *0.  129.  3730.  2611.  -0.30 

26  5 2696  1683.  *.  *0.  125.  3730.  250*.  -0.33 

28  * 2997  1578.  *.  *0.  119.  2980.  21**.  -0.28 

28  3 2398  1507.  5.  *u.  113.  2980.  2126.  -0.29 

23  2 2399  1*79.  5.  *0.  106.  2980.  2133.  -P.2P 

28  1 2900  1*95.  5.  *0.  99.  29o0.  1865.  -0.37 

26  -2  2902  1651.  5.  *0.  85.  2550.  16*3.  -0.36 

26  -3  2903  1780.  5.  *0.  79.  2550.  I960.  -0.22 

23  -*  290*  193*.  *.  *0.  7*.  2550.  1631.  -0.36 

30  1*  3093  3229.  1.  *0.  *2.  3730.  2*06.  -0.35 

30  13  309*  2932.  2.  *0.  *?.  3730.  2320.  -0.38 

30  12  3095  2t 37.  2.  *0.  -2.  3730.  2370.  -0.36 

30  11  3096  2*93.  2.  *0.  *3.  3730.  2310.  -0.38 

3C  10  3097  22*6.  2.  *0.  *3.  3730.  219*.  -0.*1 

_>0  9 3093  2012.  3.  *0.  *3.  3730.  2131  . -0.*3 

30  3 3099  1731.  3.  *2.  **.  3730.  1999.  -0.*6 

30  7 3100  156*.  3.  *0.  *6.  2930.  10*9.  -0.35 

32  5 3328  1308.  *.  *1.  **.  2980.  1796.  -0.*0 

32  * 3329  113*.  *.  *1.  *6.  2550.  1995.  -0.22 

33  9 3*03  2321.  3.  *1.  122.  2930.  1993.  -0.36 

33  3 3*0*  2086.  3.  *1.  12*.  2980.  19j0.  -0.35 

23  7 3*05  1353.  *.  *1.  127.  2980.  1733.  -0.*2 

33  6 3*06  16*2.  *.  *1.  131.  2930.  1335.  -0.38 

33  5 3*07  1**3.  5.  *1.  135.  3730.  21**.  -0.*3 

33  * 3*09  1271.  5.  *1.  1*0.  3730.  2060.  -0.*5 

33  3 3*09  1131.  6.  *1.  1**.  3730.  2087.  -0.** 

33  1 3*11  1017.  7.  *1.  1*2.  3730.  2201.  -0.*1 

35  3 3619  558.  5.  19.  *9.  2190.  1562.  -0.29 

35  2 3620  *17.  7.  19.  31.  1890.  1723.  -0.09 

**  7 **91  1582.  6.  1 5.  169.  *500.  3337.  -0.26 

**  6 **82  1333.  7.  15.  170.  *500.  3109.  -0.31 

**  5 **83  1035.  9.  15.  171.  *500.  3071.  -0.32 

**  * **  9*  9*3.  12.  15.  171.  5260.  30*9.  -0.*2 

**  3 **95  603.  18.  15.  170.  5260.  2 6°7.  -0.*9 

*5  9 *568  19CC.  *.  1*.  169.  *500.  2612.  -C.*2 

*5  8 *569  1638.  5.  1*.  168.  5260.  22fc9.  -0.38 

*5  7 *570  1*8*.  6.  1*.  1 6o • *500.  3255.  -0.28 

*5  6 *571  1280.  7.  1*.  Is*.  *500.  2330.  -0.37 

*5  5 *57?  1092.  9.  I*.  159.*  6260.  26*5.  -0.*9 

*5  * *673  916.  11.  1*.  152.  5260.  2622.  -C.*5 
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Appendix  C 

KOSCHMIEDER'S  EQUATION 


Koschmieder"s  equation  is  derived  here  from  a more  general  form  of 
the  transport  equation.  The  reason  for  including  a derivation  is  that  the 
assumptions  necessary  for  validity  and  the  expression  for  "equilibrium  luminance" 
for  the  atmosphere  both  occur  naturally. 

Ignoring  polarization,  the  propagation  of  brightness  (luminance)  in 
an  absorbing  and  incoherent  scattering  atmosphere  is  in  accord  with  the  well- 

g 

known  expression  (e.g.,  Middleton  ): 

/ f , , r(x')  \ -T(x) 

a = ( Bo  + J S (x  ) e dx'  ) e (C-l) 


Here,  B0  is  inherent  brightness  at  X = 0,  and  3 is  apparent  brightness  at  a 
distance  X.  S (. r)  is  a veiling  brightness  source  term  given  by: 


S = k 


(C-2) 


where  is  the  surround  brightness  for  a given  point  on  the  path,  and  A is  the 
volume  scattering  coefficient.  The  parameter  T (optical  thickness)  is  defined 
as  follows: 

x 

T (X)  = 3 (x‘)  d.  z'  (C-3) 

o 


where  /3  is  the  attenuation  coefficient.  Finally,  one  further  term,  the 

A 

scattering  phase  function, k , is  defined: 

k = k//3  (C-4) 


8.  W.E.K.  Middleton,  "Vision  Through  the  Atmosphere,"  University  of  Toronto 
Press,  1968. 


39 


Since  both  the  volume  scattering  coefficient  and  the  attenuation  coefficient 
are  directly  proportional  to  atmospheric  concentration,  k depends  on  the 
physical  nature  of  the  scattering  medium,  but  not  on  concentration. 

The  approximation  that  leads  to  Koschmieder's  equation  is  to  take 
both  k and  as  independent  of  path  position.  With  manipulation.  Equation 
(C-l)  becomes: 


a ’ B0eV(X)  + - e'TCX)) 


(C-5) 


For  an  infinitely  thick  path  (T-*—),  the  "equilibrium  brightness",^,  can 
be  defined. 


?oo  - <D  k 3S  cCjTL 


(C-6) 


and,  therefore: 


a popular  form  of  Koschmieder's  equation. 
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